Anderson transition in three and four effective dimensions 
for the frequency modulated kicked rotator 
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The critical exponents for the Anderson transition in three and four effective dimensions are 
discussed on the basis of previous data obtained for the frequency modulated kicked rotator. Without 
appeal to a scaling function they are shown to be in a satisfactory agreement with the theoretical 
relation known for them. 
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Recent experiments of the Garreau group with kicked 
cold atoms [1-3] renewed interest to the frequency modu- 
lated kicked rotator (FMKR), which had been introduced 
and studied some time ago [4-7]. The evolution of the 
quantum system is described by the unitary propagation 
operator 



Si = exp{-iHo{n)) ex.p{-iV{9,t)) 



(1) 



with quasiperiodic frequency modulation of kick poten- 
tial V{9,t) = ^(6*, 6*1, 6*2) with 6*1,2 = ^1,2* at d = 3 and 
V{9,t) = ¥{9,01,62,93) with 6ii,2,3 = wi,2,3i at d = 4. 
Here, the notations following those of [6]. 

The model with arctangent kick rotator (AKR) poten- 
tial was considered in [5, 7], while the frequency modu- 
lated kicked rotator (FMKR) corresponding to the quan- 
tum Chirikov standard map with modulated kick ampli- 
tude had been analyzed in [4, 6]. 

In view of this interest 1 reconsider in more detail the 
results presented in [6] for FMKR in effective dimensions 
d = 3,4 (data of Figs. 1,9 in [6] respectively). In the 
case d = 3 we have V{9,t) = fccos6'(l -I- ecosuJitcosLd2t) 
with fixed e = 0.75, irrational frequencies a;i.2 and ran- 
dom but fixed in time rotational phases Ho{n) (see [6] 
for detailed notations). The similar choice of V{9,t) = 
kcos9{l + e cos Ldit cos uJ2t cos Ldst), e = 0.9 is done for 
the case of d = 4 (data of Fig. 9 in [6]). Let me note 
that in the experiment [1-3] the rotational phases corre- 
spond to a free propagation with Ho{n) — with 
the classical chaos parameter K = kT. However, even for 
one modulation frequency the chaos border for destruc- 
tion of two-frequency invariant torus is very low so that 
the quantum chaotic dynamics mimics rather well ran- 
dom quantum phases H(){n). Indeed, according to Fig. 5 
in [8] the invariant torus with two spiral mean frequen- 
cies is destroyed a.t K 0.3 for e w 0.75 that is signifi- 
cantly smaller than the experimental values with K 6, 
T « 2.89, fc sa 2. On the basis of these arguments it is 
natural to expect that the FMKR with random phases 
has transition approximately at the same parameters as 
for the FMKR model with quadratic rotational phases 
(e.g. with r = 2). This was confirmed by the first nu- 
merical simulations for the case of quadratic rotational 
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FIG. 1: (Color online) Power law dependence for the inverse 
localization length 7 = I//1 and diffusion rate D obtained at 
asymptotically large times t at FMKR with e = 0.75 (in a 
vicinity of critical point t > 10®). Data are taken from Fig.l 
of [6] and are plotted for the fixed value of critical parameter 
kc — 1.8 for dimension d = 3 with full circles for 7 and open 
circles for D (dotted lines are drown to adapt an eye). The 
straight lines show the power law fits with fixed fee and 1/ = 
1.537 ± 0.0539, log^g 70 = -0.444 ± 0.0192 and s = 1.583 ± 
0.0511, logio Do = 0.175 ± 0.0301. 



phases performed in the proposal of Garreau experiment 
in 2005 [9] (done for effective d = 3 at e = 0.75, T = 2 
with estimated critical kc ~ 1.8). Hence, the both mod- 
els are rather similar. The transition border found in nu- 
merical simulations [1-3] is also in agreement with this 
statement. Indeed, according to the data of Fig.l in [1] 
one finds k = K/T = K/h = 1.88 ± 0.035 for e = 0.75 
and T = ?i = 2.89 that agrees with the critical value 
kc — \.% given in [6] for the case with random rotation 
phases. 

The data of Figs. 1,9 of [6] clearly show the presence of 
Anderson transition at kc ~ 1.8 (d — 3) and kc ~ 1.15 
(d = 4). These data give the inverse localization length 
7 = 1/Zi extracted from a steady-state distribution at 
asymptotically large times t > 1/7'' in the localized phase 
k < kc- In the same way in the metallic phase they give 
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FIG. 2: (Color online) Power law dependence for the inverse 
localization length 7 = l/h and diffusion rate D obtained 
at asymptotically large times t at FMKR with e = 0.9 (in a 
vicinity of critical point t > 10®). Data are taken from Fig. 9 
of [6] and are plotted for the fixed value of critical parameter 
kc = 1.15 for dimension d = 4 with full circles for 7 and 
open circles for D (dotted lines are drown to adapt an eye). 
The straight lines show the power law fits with fixed kc and 
y = 1.017±0.041, logioTo = -0.100±0.0266 and s = 2.003± 
0.0740, logio Do = -0.324 ± 0.0317 (the lowest value of D is 
not taken into account in the fit since the evolution time was 
smaller than 1/_D*). 



±12%. At the same time the statistical accuracy remains 
approximately on the same level. The fit of data by 3- 
parameter power law gives: Do = 1.626 ± 0.086, kc = 
1.728 ±0.039, s = 2.074±0.177and7o = 0.135 ± 0.0269, 
kc = 2.220 ± 0.0783, = 2.626 ± 0.145 for d = 3. 
Respectively, such a fit for d = 4 data gives Dq = 
0.397 ± 0.090, kc = 1.033 ± 0.082, s = 2.697 ± 0.263 and 
70 = 0.821±0.028, kc = 1.180±0.033, u = 1.176±0.091. 
These results show that the exact computation of the 
critical exponents s and i' remains very hard task even 
if the FMKR model is much more efficient comparing 
to the transfer matrix techniques and other numerical 
methods. It is possible that the scaling methods used in 
[1-3] correspond to a better accuracy of the exponents. 
However, this approach uses extrapolation methods com- 
bined with the scaling function which are not used in the 
data presented here. 

A separate note should be done for the AKR model. 
The results presented in [7] show very strong deviation 
from the scaling relation between s and v for 4 < d < li- 
lt is possible that at large d the critical kick parameter 
becomes relatively small kc ^ 1/d and thus this model 
becomes too close to a model with practically decoupled 
quasi-periodic driving in time that can make it rather 
specific or can require to study very small vicinity of kc 
for extraction of correct critical exponents. 



the diffusion rate D computed on asymptotically large 
times ta > from a gaussian distribution over sites. 

The results are averaged over 100 disorder realizations 
for k < kc and 10 realizations for k > kc with t ~ 10^. In 
contrast to the approach used in [1-3], the extraction of 
such asymptotic values used in [6] does not rely on the 
existence of scaling function. 

The data of [6] can be used for extraction of the scaling 
exponents v, s in the localized 7 = 7o|A:c— fc|'^ and metallic 
D = Do\kc — k\^ phases. For d — 3 the fit of data for 7 
with fixed critical point kc — 1.8 is shown in Fig. 1. In 
a similar way the fit for d = 4 is shown in Fig. 2. The 
obtained values of v and s are in a satisfactory agreement 
with the scaling relation s = (d— 2)i/ (see. e.g. [10]) both 
for d = 3, 4. For d — 3 the values of s, ^ are compatible 
with those obtained in [1-3]. According to [11] D.Delande 
obtains for d — A the critical exponents compatible with 
the scaling relation and values similar to those of Fig. 2. 

Even if formally the statistical errors are relatively 
small the actual values of v and s remain rather sen- 
sitive to variation of kc- Thus for d = 3 a variation of 
kc by ±0.05 gives variation of s by ±0.35, —0.40 and of 
V by —0.17, ±0.16. In a similar way for d = 4 a varia- 
tion of kc by ±2% gives variation of s by =f6% and v by 
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